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14.  ABSTRACT 

Photons  of  the  appropriate  energy  can  induce  photofission  in  fissile  material.  We  are  investigating  the  applicability  of  this  mechanism,  using 
photons  from  bremsstrahlung,  for  long-range  (~  100  m)  detection  of  fissile  material  using  high-power  (~  1  TW),  pulsed  technology-so  called 
intense,  pulsed  active  detection.  In  this  report,  we  describe  the  results  of  laboratory,  proof-of-principle  experiments,  supported  by  computations 
and  analyses,  in  which  a  single  pulse  of  8-MeV  endpoint  bremsstrahlung  from  the  Naval  Research  Laboratory  Mercury  pulsed-power  generator 
(200-kA  peak  current,  50-ns  pulse  width)  induces  photofission  in  depleted  uranium.  Fission  products  are  measured  using  He-3  proportional 
counters,  and  plastic  (BC408),  sodium-iodide  (NaI:Tl),  and  bismuth-genninate-oxide  scintillators.  Prompt  neutrons,  delayed  neutrons,  and 
delayed  gamma-rays  from  the  induced  photofission  are  measured  unambiguously.  Delayed  neutron  and  gamma-ray  measurements  are  also 
carried  out  with  various  thicknesses  of  low-  and  high-atomic-number  materials  attenuating  the  bremsstrahlung  and  fission  signatures.  Results 
of  simulations  and  analyses  carried  out  in  support  of  the  experiments  are  in  agreement  with  the  measurements.  In  addition,  we  demonstrate  that 
relatively  simple  variations  of  the  diode  geometry  can  increase  the  efficiency  of  inducing  fissions  three-fold  over  the  original  geometry.  These 
geometry  changes  cause  the  electron  beam  to  be  more  normally  incident  onto  the  tantalum  converter  and  thereby  increase  the  forward-directed 
bremsstrahlung  intensity  without  changing  the  driving,  electrical  power  pulse.  The  results  of  this  work  and  future  directions  for  this  research 
are  presented  in  this  report. 
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1.  INTRODUCTION 


Passive  and  active  detection  of  hidden  fissile  material  is  a  subject  of  intense  interest.1,2  In 
active  detection,  the  strength  of  the  characteristic  radiations  associated  with  naturally  occurring 
fission  are  significantly  enhanced  so  that  these  signatures  are  more  readily  measured  compared 
with  the  passive  approach  that  measures  what  can  be  relatively  weak  and  easily-attenuated 
radiations  associated  with  the  natural  decay  (i.e.,  a-decay)  of  fissile  material.2,1'4  One  approach  to 
active  detection  utilizes  the  bremsstrahlung  produced  when  an  electron  beam  strikes  a  high- 
atomic-number  material.  Photons  of  about  6  MeV  and  higher  can  induce  photofission  in  fissile 
material  to  produce  detectable  fission  products.  This  approach  was  pioneered  ’  ’  ’  ’  ’  and  is  still 
being  pursued10  using  a  pulse  train  of  bremsstrahlung  pulses  produced  by  a  linear  accelerator 
(linac),  commonly  used  for  medical  applications.  The  linac  approach  is  attractive  because  it 
employs  nearly-commercial-off-the-shelf,  compact  technology.  However,  because  of  the 
relatively  low  current  associated  with  linacs  (usually  less  than  1  mA  average  current),  it  can  take 
several  minutes  of  bremsstrahlung  irradiation  at  relatively  high  electron  endpoint  energy  (20  MeV 
or  higher)  to  induce  sufficient  fissions  for  long -standoff  detection  applications  (~  100  m).  The 
relatively  long  measurement  time  can  be  a  disadvantage  in  a  practical  application.  While  the  high 
endpoint-energy  of  the  linac  approach  results  in  more  photons  emitted  per  electron  near  the 
maximum  of  the  photofission  cross  section  (~  14  MeV)11  than  with  moderate-endpoint-energy 
photons  (e.g.,  8-12  MeV),  the  resulting  high-energy  photons  can  induce  collateral  photo¬ 
activation  of  the  environment  and  other  benign  materials  in  the  photon  beam  path.12,13  The 
resulting  induced  (active)  background  radiation  can  severely  complicate  a  high-confidence 
measurement  of  the  fission  signature.  Also,  a  significant  natural  (passive)  background  (e.g.  from 
cosmic  rays)  can  accumulate  with  long  measurement  time. 

An  alternative  to  the  linac  approach  relies  on  the  use  of  TW-level  pulsed-power  technology, 
similar  to  that  used  in  laboratory  simulation  of  nuclear  weapon  effects.14  The  complete 
measurement  of  a  fission  signature  made  after  a  single  pulse  is  generally  termed  intense,  pulsed, 
active  detection  (IPAD).15  The  high-power  electrical  pulse  is  used  to  create  a  short  (~  100  ns), 
intense,  bremsstrahlung  pulse  that  induces  photo-fission  in  fissile  material.15, 16,17,ls  The  IPAD 
approach  has  the  potential  advantage  of  being  able  to  measure  within  a  few  seconds  after  the 
pulse  both  the  prompt  neutron  signature  as  well  as  the  beta-delayed  neutron  and  y-ray  signatures 
associated  with  the  induced  photo-fissions.  Prompt  neutrons  may  be  detected  within  a  few 
hundred  nanoseconds  after  the  irradiation  pulse,  depending  on  the  detector  distance.  The  much 
shorter  measurement  time  of  the  IPAD  technique  reduces  the  natural  background  signal  to  a  much 
smaller  fraction  of  the  measured  signal  than  obtained  with  a  conventional,  repetitively-pulsed 
linac  approach.15  However,  progress  is  required:  (1)  to  develop  detectors  that  recover  promptly 
from  the  scattered  radiation  associated  with  the  interrogating  pulse,  (2)  to  minimize  the  induced 
background  (i.e.,  by  using  moderate-energy  bremsstrahlung),12’13’15  and  (3)  to  make  pulsed-power 
driver  technology  more  robust  and  compact.19,20 

Earlier  IPAD  experiments21,22  used  6-  to  7-MeV  characteristic  gamma-rays  produced  via  the 
19F(p,ay)160  nuclear  reaction  to  induce  fission  in  a  depleted  uranium  (DU)  object.23  An  intense 
characteristic  y-ray  pulse  was  produced  with  the  Naval  Research  Laboratory  Mercury  inductive 
voltage  adder  (IV A)24  operated  in  positive  polarity25  to  generate  an  intense  proton  beam26  that 
bombarded  a  fluorinated  target  (CF2).  The  characteristic  gamma-rays  induced  fission  in  a  DU 
object,  and  fission  neutrons  were  measured.  The  results  are  reported  in  Refs.  21  and  22.  While 
products  from  induced  fission  were  measured  using  this  technique,  the  low  efficiency  of 
characteristic  gamma  ray  production  and  their  inherent  isotropic  emission  preclude  application  of 
this  approach  to  long  distance  (~  1 00-m)  standoff. 

Early  experimental  investigations  of  the  IPAD  approach  with  moderate  energy 
bremsstrahlung  and  analyses  of  these  experiments  are  presented  in  this  report.  For  these 
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experiments.  Mercury  was  modified  to  increase  its  operating  voltage  from  6  to  8  MV  (with  200- 
kA  peak  current  and  50-ns  electrical  pulse  width).27  The  8-MeV  bremsstrahlung  from  Mercury 
irradiated  a  30-cm  x  30-cm  x  2.5-cm-thick  DU  plate  from  which  fission  products,  including 
prompt  neutrons,  delayed  neutrons,  and  delayed  y-rays,  where  measured.16'17'18  Delayed  y-rays 
were  measured  with  bismuth-germinate-oxide  (BGO)  and  sodium-iodide  (NaI:Tl)  scintillator 
detectors.  Measurement  of  a  prompt-y  fission  signature  was  assumed  unachievable  because  of 
the  much  more  intense  scattered  bremsstrahlung  in  the  same  spectral  band  and  essentially  at  the 
same  time.  Delayed  neutrons  were  measured  with  He-3  proportional  counters.  Prompt  neutrons 
were  measured  with  a  fast-time -response  plastic-scintillator  detector  to  separate  neutrons  from 
photons  (scattered  bremsstrahlung)  by  virtue  of  their  time-of-flight  difference.  Unambiguous 
signatures  of  induced  fission  were  obtained  from  the  delayed  y-rays  measurements  and  from  both 
the  delayed  and  prompt  neutron  measurements. 

Analytic  and  computational  analyses  were  used  to  interpret  the  measurements  and  guide  the 
experiments.  The  large-scale -plasma  (LSP)  particle -in-cell  code28  was  used  to  model  the 
dynamics  of  electron  flow  in  the  Mercury  IVA  and  the  electron  beam  (e-beam)  dynamics  in  the 
diode.  The  integrated  tiger  series  (ITS)  Monte  Carlo  electron-photon  transport  code29  was  used  to 
model  bremsstrahlung  production  and  transport.  The  Monte-Carlo  N-particle  extended 
(MCNPX)  code30  was  used  to  model  relevant  nuclear  reactions,  neutron  generation,  and  neutron 
transport.  The  measurements  are  also  useful  in  benchmarking  these  codes.16  18 

Measurements  of  delayed  neutrons  and  y-rays  from  induced  fission  also  were  carried  out  with 
materials  interspersed  between  the  DU  plate  and  the  bremsstrahlung  source.  Various  thicknesses 
of  low-  and  high-atomic-number  (low-  and  high-Z)  materials  were  used  in  these  studies,  as  well 
as  composites  of  low-  and  high-Z  materials.  Fission  signals  equal  to  or  greater  than  3  times  the 
background  standard  deviation  are  obtained  with  up  to  160  g/cm2  areal  mass  density  of  high-Z 
attenuating  material.  Simulations  of  the  fission  signatures  after  radiation  transport  through 
various  areal  mass  densities  of  attenuating  material  are  in  general  agreement  with  the  measured 
fission  signatures. 

In  addition  to  these  fission  signature  measurements,  other  experiments  demonstrate  an 
increase  in  the  bremsstrahlung  intensity  on  the  DU  plate  by  a  factor-of-three.  This  is 
accomplished  by  adjusting  the  anode-cathode  (AK)  gap  and  the  diameter  of  the  outer  conductor 
in  the  e-beam-diode  region  to  make  the  electron  trajectories  striking  the  tantalum  converter  more 
normally  incident,  thereby  increasing  the  radiation  intensity  in  the  forward  direction. 

All  of  these  results  are  discussed  in  more  detail  in  the  following  sections  of  this  report. 
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2.  RADIATION  SOURCE  AND  MERCURY  PULSED  POWER 


The  design  and  modeling  of  the  bremsstrahlung  diode  are  described  in  detail  in  Ref.  17.  The 
front  end  of  the  Mercury  IVA  is  illustrated  in  Fig.  1 ,  where  the  outputs  of  the  generator  cavities 
(not  shown)  are  connected  to  the  load  via  a  magnetically  insulated  transmission  line  (MITL).  ’1  A 
brief  description  of  the  initial  diode  geometry  follows.  The  17-cm  diameter  cylindrical  MITL 
center  conductor  on  the  Mercury  generator  is  at  negative  high  voltage  and  is  terminated  in  a  17- 
cm  diameter  hemisphere.  A  relatively  large  axial  AK  gap  is  used  to  minimize  electron  emission 
from  the  hemisphere  so  that  the  load  presents  a  high  impedance  to  the  machine  and  the  MITL 
fixes  the  effective  load  impedance  to  its  self-limited  impedance.  ’  ’  Thus,  the  voltage  is 
fixed  at  8  MV  and  a  significant  fraction  of  the  electron  current  incident  on  the  anode  tantalum 
converter  comprises  MITL  flow  electrons  (or  vacuum-flow  current)  emitted  from  along  the  MITL 
center  conductor  in  what  is  called  “self-limited”  or  “line-limited”  flow.34'35  In  Fig.  1,  the 
38-cm-diameter  converter  is  located  23  cm  from  the  tip  of  the  hemisphere  on  the  center 
conductor.  Because  of  this  large  axial  gap,  the  end  of  the  MITL  acts  as  a  radial  diode  emitting 

the  remainder  of  the  load 
current.  This  diode-emitted 
electron  current  joins  the 
MITL  flow  current  and  this 
combined  total  electron 
current  strikes  the  converter. 
All  the  data  described  in 
Secs.  3  through  5  were  taken 
with  this  23-cm  diode  AK 
gap.  The  data  in  Sec.  6  were 
taken  with  a  40-cm  AK  gap. 
Section  VII  describes  diode 
improvement  work  where  the 
AK  gap  is  varied  to  optimize 
the  electron  angles  of 
incidence  on  the  Ta 
converter.  The  Ta  converter 
thickness  in  all  of  the 
experiments  is  only  1 50- pm 
(much  less  than  optimum  for 
maximum  dose)  and  has  an 

aluminum  beam-stop.  This  converter  geometry  optimizes  the  ratio  of  higher  energy  photons  that 
induce  photofission  to  lower  energy  photons  that  produce  unwanted  dose.17 

To  minimize  the  possibility  of  ion-beam  production  in  the  diode,  the  large -radius  hemisphere 
is  used  to  terminate  the  MITL  so  as  to  reduce  the  emitted  e-beam  current  density  and  associated 
deposition  heating  of  the  anode  surface.  The  converter  is  also  covered  with  a  thin  (18-pm  thick) 
bare  A1  front  layer  both  because  A1  has  a  higher  thermal  threshold  for  anode  plasma  formation 
and  little  beam  energy  is  deposited  within  its  thin  layer.  By  avoiding  anode  plasma  formation, 
which  would  result  in  ion-beam  generation  and  partial  space-charge  neutralization  of  the  e-beam, 
minimal  e-beam  focusing  (or  pinching)  occurs  in  the  diode,  as  desired.  This  front  A1  layer  also 
prevents  ions  that  form  on  the  Ta  converter  from  entering  the  vacuum  diode  region.  However,  as 
discussed  below,  these  procedures,  while  preventing  anode  plasma  formation,  do  not  totally 
eliminate  ions  from  the  diode  region. 

To  provide  multiple  reproducible  shots  without  opening  the  system,  the  cathode  and  outer 
anode  conductors  were  not  coated  with  the  usual  thin  graphite  layer  (e.g.,  using  Aerodag36),  but 
instead  are  bare  A1  surfaces.  Multiple  shots  with  reproducible  current  and  voltage  waveforms  are 


Figure  1.  Schematic  of  the  Mercury  IVA  front  end  and 
bremsstrahlung  diode.  The  AK  gap  shown  is  23  cm. 
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achieved  without  opening  the  system  as  reported  in  Ref.  17.  A  thick  steel  collimator  (20-cm  thick 
upstream  of  the  converter,  30-cm  thick  downstream  of  the  converter)  surrounds  the  diode  to 
attenuate  off-axis  bremsstrahlung  (see  Fig.  1).  This  collimator  has  a  33-cm-diameter  downstream 
bore  and  reduces  the  off-axis  dose  (angles  >  20°)  by  more  than  an  order  of  magnitude.17 

Diode  performance  including  load  current,  load  voltage,  and  bremsstrahlung  production  are 
described  in  Ref.  17.  The  electron  dynamics  of  the  diode  region,  including  both  vacuum-flow 
electrons  from  along  the  length  of  the  MITL  and  electrons  generated  in  the  radial  diode  at  the  end 
of  the  MITL  (see  above),  are  modeled  with  the  LSP  particle -in-cell  code,  and  the  results 
compared  with  various  measurements.  As  shown  in  Fig.  2,  the  calculated  (total)  load-current, 
I io;i(i,  and  load-voltage,  Vioad,  waveforms  are  in  agreement  with  the  measured  waveforms.  These 
quantities  are  measured  at  the  end  of  the  MITL  before  the  hemisphere.  Also  shown  in  Fig.  2  is  a 
measured  bremsstrahlung  pulse  associated  with  the  electrical  pulse.  In  the  simulations,  the 
electron  flow  in  the  diode  is  modified  by  a  small  amount  of  ion  charge  created  at  the  anode 

surface  by  stimulated  emission 
resulting  from  electron  impact.’7 
This  stimulated  ion  emission  in  a 
given  area  is  determined  by 
setting  an  ion  emission  fraction 
in  terms  of  ions  emitted  per 
incident  electrons  in  that  area. 
Stimulated  ion  emission  fractions 
ranging  from  0.2%  to  0.6%  are 
shown  to  have  a  substantial 
impact  on  the  electron 
distribution  and  angles  of 
incidence  at  the  converter,18  and, 
as  a  result,  on  the  forward-going 
bremsstrahlung  intensity.  The 
issue  of  the  ion  fraction  will  be 
discussed  further  in  Sec.  4, 
below. 

The  dose  associated  with  the 
bremsstrahlung  pulse  shown  in 
Fig.  2  (and  for  the  23-cm  AK 
gap)  at  5  m  from  the 
bremsstrahlung  converter  is 
typically  20  rad  (CaF2).  Various 
simulations  (discussed  in  what  flows)  suggest  that  this  bremsstrahlung  pulse  induces  about 
2.25xl09/[RSo(m)]2  fissions  in  the  DU  plate  when  it  is  a  distance  Rso(m)  from  the  bremsstrahlung 
converter. 
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Figure  2.  Comparison  of  experimental  (thick  lines)  and  LSP- 
computed  (thin  lines)  load  voltage  and  total  current  waveforms 
for  the  Mercury  bremsstrahlung  diode.  Also  shown  is  the 
experimental  photodiode  x-ray  waveform  (arbitrary  scale). 
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3.  DELAYED  y-RAY  MEASUREMENTS 
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The  arrangement  for  the  delayed-gamma-ray  measurements  is  illustrated  in  Fig.  3.  The  DU 
plate  is  located  along  the  Mercury  center  line  5  m  from  the  bremsstrahlung  source  but  tilted  at  a 
45°  angle  to  center  line.  As  shown  in  Fig.  3,  the  BGO  detector  is  located  2.65  m  from  the  DU 
plate  and  consists  of  six  5-cm  diam  by  5-cm  thick  BGO  scintillators  with  photomultipliers 
mounted  in  a  3  x  2  array.  Each  detector  was  collimated  with  lead  and  30%-borated  polyethylene 
to  partially  mitigate  stray  photon  and  room-scattered  neutron  background.  Individual  pulses  are 

recorded  in  list  mode.  The 
bremsstrahlung  _  pulses  can  be  summed  in 

various  time  bins  to 
provide  the  count  rate  as  a 
function  of  time.  Also,  the 
pulse  height  is  used  to 
separate  the  counts  into 
two  broad  y-ray  energy 
ranges:  a  broad,  low- 
energy-cutoff  range  with 
Ey  >  350  keV,  and  a  high- 
energy  range  with  3  MeV 
<  Ey  <  6  MeV.  The  350- 
keV  cutoff  energy  is 
selected  to  avoid  detecting 
the  140-keV  y-ray  from 
75mGe  (47-s  half-life), 
which  is  produced  by 
neutron  activation  in  the 
BGO  detector.  The  3- 

MeV  cutoff  energy  is  selected  to  eliminate  the  naturally  occurring  2.614-MeV  y-ray  from  2()8T1  (in 
the  decay  chain  of  208Th). 

Comparisons  of  measurements  with  the  DU  plate  and  a  lead  surrogate  are  presented  in  Fig.  4 
for  the  broad,  low-energy-cutoff  y-rays  (Figs.  4a  and  4c)  and  for  the  high-energy  y-rays  (Figs.  4b 
and  4d).  In  Figs.  4a  and  4b,  the  counts  in  each  energy  range  are  summed  in  2-ms  bins  from  1  ms 
to  0. 1  s  and  are  plotted  as  a  function  of  time  at  the  center  of  each  bin.  In  Figs  4c  and  4d,  the 
counts  in  each  energy  range  are  summed  in  20-ms  bins  from  1  ms  to  1  s  and  are  plotted  as  a 
function  of  time  at  the  center  of  each  bin.  The  error  bars  on  the  counts  in  these  figures  are 
statistical  uncertainties.  At  early  time  in  Figs.  4a  and  4b,  both  the  DU  and  the  lead  counts 
increase  as  the  high  voltage  on  the  photomultiplier  tube  begins  to  recover  from  the  intense 
bremsstrahlung  pulse.  After  this  4-  to  6-ms  detector  recovery,  both  the  low-energy-cut  and  the 
high-energy-range  counts  decrease  in  similar  fashion.  At  this  time,  however,  the  high-energy- 
range  counts  are  about  four-to-five  times  smaller  than  the  low  energy  counts.  This  is  expected 
because  of  fewer  photons  and  lower  detector  efficiency  for  higher-energy  delayed  y-rays. 
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Figure  3.  Top  view  of  the  bremsstrahlung  diode,  DU  plate  and  BGO 
detectors  for  delayed  y-ray  measurements. 
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As  seen  in  Figs.  4(a)  and  (c),  the  induced  (active)  background  measured  with  the  lead 
surrogate  in  the  first  20  to  30  ms  masks  the  beta-delayed  y-ray  signature  from  the  DU  fission 
products  for  both  the  broad-low-energy-cutoff  and  high-energy-range  cases.  This  background 
persists  for  too  long  to  be  associated  with  room-scattered  bremsstrahlung  from  the 
bremsstrahlung  interrogating  pulse  (~  30-ns  FWHM,  see  Fig.  2).  However,  bremsstrahlung  can 
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Figure  4.  Delayed  g-ray  counts  plotted  for  the  low  energy  cutoff  in  (a)  and  (c)  and  for  the  high-energy 
range  in  (b)  and  (d).  The  red  squares  are  for  shot  821  with  the  DU  plate,  and  the  black  circles  are  for  shot 
823  with  the  lead  surrogate.  The  counts  are  summed  in  2-ms  time  bins  in  (a)  and  (b)  and  in  20-ms  time 
bins  in  (c)  and  fd). 


interact  with  the  local  environment  to  produce  neutrons  through  photonuclear,  (y,n),  reactions. 
These  neutrons  may  then  thermalize  and  be  captured  by  nuclei  in  the  environment.  These  nuclei 
subsequently  undergo  (n,y)  reactions  that  produce  short-lived  y-ray  activity  in  materials 
comprising  the  entire  environment  as  well  as  material  surrounding  the  detector.  Also,  about  I  ()  ’ 
neutrons,  produced  in  the  e-beam  diode,  may  contribute  to  this  neutron-induced,  background  y- 
ray  activity.17  This  process  can  last  for  as  long  as  the  photo  neutrons  take  to  thermalize,  ~  10  ms 
(or  longer)  similar  to  the  time  scale  indicated  in  Figs.  4(a)  and  (c)  for  the  active  delayed-gamma 
background. 

The  DU  and  lead  surrogate  measurements  are  compared  at  later  time  using  the  20-ms  time 
bins  in  Figs.  4(c)  and  (d).  After  about  40  ms,  the  DU  counts  plateau  at  about  200  counts  per  bin 
for  the  low-energy-cutoff  case  (Fig.  4(c))  and  less  than  10  counts  per  bin  for  the  high-energy- 
range  case  (Fig.  4(d)).  The  sums  of  the  low-energy-cutoff  counts  from  40  ms  to  1  s  are  7840 
counts  with  DU  and  49  counts  with  lead,  corresponding  to  a  ratio  of  160.  For  the  high-energy 
cut,  these  sums  over  the  same  time  period  are  280  counts  with  DU  and  no  counts  with  lead. 
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Clearly,  unambiguous  detection  of  the  DU  can  be  achieved  within  1  s  using  a  single,  intense 
pulse.  The  low  background  associated  with  the  high-energy  y-rays  can  be  advantageous  when 
quantifying  detection  confidence. 
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4.  DELAYED  NEUTRON  MEASUREMENTS 


Delayed  neutrons  are  detected  with  an  array  of  five  He-3  proportional  counters.  The  setup 
for  these  measurements  is  illustrated  in  Fig.  5.  These  detectors  are  encased  in  polyethylene, 
cadmium,  and  borated  rubber  to  emphasize  their  response  to  the  delayed  neutron  fission  spectrum 

relative  to  lower 
energy  neutrons.  The 
design  and  shielding 
of  these  detectors  are 
described  in  detail  in 
Refs.  22,  38,  and  39. 
As  shown  in  Fig.  5, 
the  five  detectors  are 
located  at  distances 
ranging  from  1  m  to 
4.5  m  from  the  DU 
plate  and  the  DU  plate 
is  on  the  Mercury 
center  line  5  m  from 
the  bremsstrahlung 
source  and  tilted  at  a 
45°  angle  to  the  center 
line.  Large  pulses  are 
produced  by  low 
energy  neutrons  due  to 
the  positive  Q-value  (0.762  MeV)  of  the  3He(n,p)3H  reaction  in  the  detector. 

The  first  100  ms  of  neutron  pulses  after  the  bremsstrahlung  pulse,  measured  using  the  He-3 
detector  located  1  m  from  the  front  face  of  the  DU,  are  compared  in  Fig.  6  for  shots  with  DU, 
with  lead  only,  and  with  no  object.  For  these  data,  the  detectors  recover  from  the  bremsstrahlung 
flash  in  10  ms  or  less.  Many  neutron  pulses  are  measured  with  the  DU,  but  no  neutron  pulses  are 
measured  with  either  the  lead  surrogate  or  no  object.  For  this  case,  summing  the  pulses  recorded 
with  the  DU  in  100  ms  could  provide  a  clear  delay ed-neutron  signal  with  essentially  no  active 
background.  This  time  scale  is  so  short  that  the  passive  background  is  negligible  (the  passive 
background  count  rate  at  Mercury  for  this  type  of  detector  is  about  0.024  counts/s).  This  low 
passive  background  contribution  to  the  measurement  is  an  inherent  advantage  of  the  IPAD 
approach. 

Figure  7  shows  a  comparison  of  the  counts  recorded  on  the  He-3  detector  located  1  m  behind 
the  DU  (see  Fig.  5)  summed  in  30-ms  time  bins  as  a  function  of  time  for  a  shot  with  the  DU  plate, 
a  shot  with  the  lead  surrogate,  and  a  normalized  fit  of  the  6-exponent  decay  model  for  DU 
photofission.40  The  data  with  lead  show  one  count  early  in  time  (possibly  from  lead  photo¬ 
neutrons)  and  one  late  in  time.  The  DU  data  are  in  agreement  with  the  decay  model. 


Figure  5.  Top  view  of  the  arrangement  of  five  He-3  detectors,  the  DU 
plate,  and  the  bremsstrahlung  diode  for  delayed  neutron  measurements. 
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Figure  7.  Signals  from  the  He-3  detector  located  1  m  from  the 
front  face  of  the  DU  plate.  Measurements  with  the  DU  plate 
(shot  811),  the  lead  surrogate  (shot  812)  and  no  object  (shot  813) 
are  shown. 


In  Fig.  8,  the  calculated 
neutron  fluence  incident  on  the 
He-3  detectors  as  a  function  of 
distance  from  the  DU  plate  (see 
Fig.  5)  is  compared  with  the 
neutron  fluence  inferred  from 
the  measured  pulses.  The 
bremsstrahlung  spectrum 

incident  on  the  DU  plate  is 
calculated  with  LSP/ITS 
assuming  a  constant  0.4%  ion 
fraction  for  stimulated  emission 
(see  Sec.  2).  Then  using  the 
photon  spectrum  as  input, 
MCNPX  is  used  to  calculate  the 
neutron  spectrum  emitted  from 
the  DU  and  transported  to  each 
He-3  detector.  This  calculation 
includes  the  delayed  neutrons 
from  the  induced  fission  in  the 
DU  plate  as  well  as  the 
contribution  from  the  (induced- 
fission)  neutrons  that  reflect  off 
the  composition-specific 

concrete  walls,  floor,  and  ceiling 


of  the  Mercury  cell.16  These 
calculations  are  represented  by  the 
black  squares  in  Fig.  8.  The  He-3 
counts,  measured  over  only  100  s, 
are  corrected  to  saturation  at  200  s 
using  the  6-exponent-decay  model40 
and  averaged  over  three  shots  for 
each  detector  distance.  MCNPX  is 
used  again  to  convert  these  measured 
He-3  counts  into  inferred  incident 
delayed-neutron  fluence.16'13  for 
comparison  with  the  computed 
fluence.  The  previously-calculated 
spectrally  dependent  incident  neutron 
fluence,  including  the  scattered 
delayed-neutron  contribution,  must 
be  used  for  this  conversion  as  the 
absolute  detector  response39  is 
neutron-energy  dependent.  These 
data  are  plotted  as  red  diamonds  in 
Fig.  8.  Folding  in  the  detector  area,22 
each  He-3  detector  count  at  1  -m  from 
the  DU  plate  is  associated  with  about 
8.9xl04  induced  fissions  in  the  DU 
plate.  The  error  bars  on  these  data 


Figure  6.  He-3  detector  counts  in  30-ms  bins  vs.  time 
for  shot  821  with  the  DU  plate  (red  squares),  and  shot 
823  with  the  lead  surrogate  (blue  circles).  The  DU  data 
are  compared  with  the  6-component-decay  model  for 
photofission.  Representative  error  bars  indicate 
statistical  uncertainty.  The  lead  data  show  a  total  of  2 
counts:  one  before  0.03  s  and  one  at  74  s  (i.e.,  most  data 
points  are  at  zero). 
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include  the  quadratic 
combination  of  the  statistical- 
counting  uncertainty  and  the 
standard  deviation  of  the 
average  of  three  shots. 

The  scaling  of  these  data  is 
not  inverse-square  with 
distance.  For  example,  the 
detector  at  4.5  m  measures 
more  than  twice  the  counts 
expected  by  inverse-square 
scaling  from  the  detector  at  1 
m.  Calculations  of  neutron 
scattering  and  absorption  by 
the  floor,  walls,  and  ceiling 
using  the  MCNPX  code 
provide  an  explanation  for  the 
measured  scaling  with 
distance.16 

Results  from  calculations 
using  two  different 
photofission  cross  section 
sources  also  are  shown  in  Fig. 
8.  MCNPX  uses  cross  sections 
from  Blohkin,41  calculated  with  a  nuclear  evaporation  model.  Calculations  using  measured  cross 
sections  from  Caldwell,11  shown  as  blue  triangles  in  Fig.  8,  are  in  closer  agreement  with  the 
measurements.  At  photon  energies  below  9  MeV,  the  Caldwell  cross  sections  are  about  two 
times  larger  than  the  Blohkin  cross  sections,  but  at  higher  energy  these  cross  sections  are  similar. 
For  the  8 -MeV  endpoint  bremsstrahlung  in  these  experiments,  the  measured  cross  sections  of 
Caldwell  give  a  better  fit  to  the  data.  An  inverse-square  dependence  supplemented  by  a  constant 
term  provides  a  reasonable  fit  to  the  results  of  measurements  and  MCNPX  simulations  (solid 
lines  in  Fig.  8).  The  inverse-square  term  represents  direct  neutron  irradiation  from  the  DU, 
whereas  the  constant  term  represents  room-scattered  neutrons.  The  parameters  for  these  fits  are 
listed  in  Table  I.  As  shown  in  Table  1,  the  neutron  fluence  from  MCNPX  is  about  a  factor  of  two 
smaller  than  the  measurements,  while  calculations  with  the  Caldwell  cross  sections  are  generally 
in  better  agreement  with  the  measurements. 

Table  I.  Fit  parameters  for  scaling  neutron  fluence  F  with  distance  R(cm). 


Data 

Source 

Equation  F  =  A  R 2  +  B 

A  (cm2) 

B(cm2) 

Measurements 

59.3  +  1.8 

5.5 +  0.9 

MCNPX 

27.8+0.8 

5.1  +0.4 

Caldwell 

38.9  +  1.0 

7.2  +  0.5 

Figure  8.  Comparison  of  the  measured  scaling  of  neutron  intensity 
as  a  function  of  distance  (red  diamonds)  with:  calculations  using 
MCNPX  with  its  usual  cross  sections  (black  squares)  and 
calculations  using  cross  sections  from  Caldwell  (blue  triangles). 
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It  should  be  reiterated  that  these  MCNPX  computations  used  a  forward-going  bremsstrahlung 
spectrum  corresponding  to  results  from  an  LSP  simulation  with  a  constant  0.4%  ion  fraction  for 
stimulated  emission.  Ion  fraction  values  of  0.2-0.4%  provide  the  best  agreement  with 
measurements  of  x-ray  dose  and  fission-neutron  yield  when  compared  with  predictions  from  ITS 
Cyltran  computations  of  x-ray  radiation  transport  to  the  DU  plate  and  detectors.18  Other  semi- 
analytic  computations  using  the  Caldwell  cross  sections  and  spectra  for  other  ion  fractions  show 
that  the  neutron  fluence  versus  distance  is  in  good  agreement  with  measurements  for  a  slightly 
smaller  ion  fraction  than  0.4%. 18  This  comparison  is  shown  in  Fig.  9  for  four  values  of  the  ion 
fraction.  The  individual  fluences  for 
the  three  shots  averaged  in  Fig.  8  are 
displayed  in  Fig.  9  with  their  values 
reduced  by  a  constant  4.5  n/cm2  to 
remove  in  an  approximate  way  the 
room-scatter  contribution. 

Results  show  that  a  delayed 
neutron  signature  is  detected  that  is 
consistent  with  simulations  and 
analysis.  The  signal  is  significantly 
above  induced  background  but  the 
induced  background  must  be  invoked 
to  explain  the  scaling  of  the  neutron 
fluence  with  DU-to-detector  distance. 

For  the  photon  energies  used  here, 
the  Caldwell  photofission  cross 
sections  give  better  agreement  with 
the  data  than  the  than  the  Blohkin 
cross  sections  used  in  MCNPX.  A 
small  ion  current  in  the  diode, 
associated  with  stimulated  ion 
emission  from  the  anode,  is  needed 
for  the  computed  neutron  fluences  to 
agree  with  the  measured  neutron 
fluences.  As  will  be  seen  in  Sec  7, 
below,  this  small  ion  current  is  also  needed  to  explain  the  angular  distribution  of  the 
bremsstrahlung-producing  electrons. 


DU  to  detector  distance  (m) 

Figure  9.  Room-subtracted  fluence  on  the  He-3  detectors 
vs  distance  from  the  DU  plate  compared  with  analysis  using 
Caldwell  cross-sections  for  0%  to  0.6%  ion  fractions.  The 
fluences  for  the  3  displayed  shots  are  averaged  and  plotted 
in  the  previous  figure 
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5.  PROMPT  NEUTRON  MEASUREMENTS 


outside  test  cell 

Figure  11.  Top  view  of  the  bremsstrahlung  diode,  the  DU  object, 
and  the  neutron  TOF  detector  for  prompt  neutron  measurements. 


Prompt  neutrons  emitted 
from  the  DU  plate  were 
measured  with  a  time-of-flight 
(TOF)  detector  mounted  within 
a  lead  shield.42  The  detector 
consisted  of  a  5.1-cm-diameter 
x  5.1-cm-thick  plastic 
scintillator  (BC408)  optically 
coupled  to  a  5.1 -cm  diameter 
photomultiplier  tube  (XP2020). 
The  arrangement  is  illustrated 
in  Fig.  10.  The  detector  is 
located  outside  the  Mercury 
test  cell  and  views  the  DU  plate 
at  90°  to  the  axis  through  a  15- 
cm-diameter  hole  in  a  2.43-m 
thick  ilmenite  concrete  wall. 
The  DU  is  located  on  the 
generator  axis  6. 17  m  from  the 
tantalum  converter  and  is 
rotated  45°  to  the  incident 
bremsstrahlung  direction  so 
that  the  neutron  detector  views 


the  front  face  of  the  plate.  The 
detector  is  located  6.17  m  from  the 
DU  plate  and  is  shielded  with  14-cm 
(front)  and  6-cm  (sides)  of  lead  to 
attenuate  the  intense  bremsstrahlung 
scattered  from  the  DU  so  that  the 
photomultiplier  tube  can  recover  in 
time  to  respond  to  the  energetic, 
prompt  neutrons  whose  arrival  to  the 
scintillator  are  delayed  by  their  TOF. 

The  DU  was  replaced  with  a  lead 
surrogate  to  compare  the  detector 
response  to  fission  neutrons  with  the 
response  to  photonuclear  neutrons 
from  (y,n)  reactions  in  the  lead 
surrogate. 

Neutron  TOF  signals  measured 
with  the  DU  plate,  the  lead 
surrogate,  and  no  object  are 
compared  in  Fig.  1 1 .  With  DU  and 
lead,  the  initial  pulses  are  of  similar 
magnitude  and  recover  to  the 
baseline  about  150  ns  after  the  rise 
(the  DU  signal  is  obscured  by  the 
lead  signal  at  this  early  time).  With 
no  object,  this  initial  pulse  is  greatly  reduced,  consistent  with  this  pulse  being  associated  with 


Figure  10.  Measured  neutron  TOF  detector  signals  for  shot 
910  with  the  DU  object  (red),  for  shot  911  with  the  lead 
surrogate  (blue),  and  for  shot  914  with  no  object  (green). 
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scattered  bremsstrahlung.  Also,  no  neutron  pulses  are  detected  without  an  object.  With  the  DU 
plate,  a  large  signal  from  piled-up  neutron  pulses  starts  about  180  ns  after  the  initial 
bremsstrahlung  rise,  corresponding  to  the  maximum  neutron  energy  from  TOF  of  about  6  MeV. 
The  initial  piled-up  neutron  pulses  are  larger  in  amplitude  than  the  later  pulses,  consistent  with 
either  higher  energy  neutrons,  or  more  coincident  neutrons  reaching  the  detector  early  in  time,  or 
both.  Fission-neutron  energies  may  extend  up  to  10  MeV,  while  photonuclear  neutrons  from  the 
238U(y,n)237U  reaction  with  a  6.2-MeV  threshold  have  energies  less  than  1.8  MeV  for  an  8-MeV 
endpoint  bremsstrahlung  spectrum.  These  measurements  confirm  that  this  TOF  detector  is 
measuring  prompt  fission  neutrons  from  the  DU.  With  the  lead  plate,  many  fewer  pulses  are 
recorded  than  with  the  DU  plate.  Most  of  these  neutrons  are  consistent  with  less  than  1.3-MeV 
photonuclear  neutrons  from  the  207Pb(y,n)  reaction,  which  has  a  6.7-MeV  threshold  energy. 

The  few  small  pulses  with  lead  between  200  and  300  ns  correspond  to  higher  TOF  energies 
and  cannot  be  so  explained,  though  they  are  observed  on  all  lead  shots.  One  speculative 
explanation  for  these  pulses  in  the  200-300  ns  range  is  neutron  TOF  from  high  energy  neutrons 
produced  in  the  diode.  From  activation  measurements  on  the  Mercury  center  conductor,  and 
rhodium  activation  foil  measurements,  we  estimate  that  approximately  1  kA  of  ions  are  generated 
by  the  high-power  e-beam  striking  the  anode  converter  package.  This  is  consistent  with  an  ion 
fraction  of  about  0.5%,  similar  to  those  discussed  in  Secs.  4  and  7.  These  ions  would 
predominantly  originate  from  hydrocarbon  contaminates  on  the  anode  surfaces.  From  isotopic 
ratios,  approximately  100  mA  of  this  ion  current  could  be  from  deuterons,  which  would  result  in 
neutrons  from  the  ”  Al(  H,n )  Si  reaction.  This  reaction  has  a  positive  Q-value  (+9.36  MeV)  and 
zero  threshold,  so  it  is  possible  for  neutrons  with  8  MeV  (diode  peak)  +  9.37  MeV  (Q-value)  = 
17.36  MeV  to  be  produced.  The  lead  surrogate  would  provide  an  excellent  scattering  location  to 
direct  some  of  these  neutrons  at  the  detector,  and  without  the  lead  surrogate  these  neutrons  would 
not  be  observed.  The  DU  plate  would  also  scatter  these  neutrons  into  the  detector  in  a  similar 
manner  to  the  lead  surrogate,  but  the  detector  response  with  DU  is  dominated  by  fission  neutrons 
at  this  time,  and  the  small  diode  neutron  signal  would  not  be  distinguishable  from  the  fission 
signature.  The  TOF  for  these  neutrons  would  be  -220  ns.  It  must  be  stressed  this  explanation  is 
a  conjecture  at  this  point  and  this  mechanism  must  be  studied  in  more  detail  to  determine  if  the 
detector  results  are  consistent  with  the  predicted  count  rate  vs  time  for  this  particular  reaction,  as 
well  as  for  any  other  neutron-producing  ion  reactions  that  would  also  occur  inside  the  diode. 

The  measured  prompt  neutron  signal  can  be  compared  with  the  calculated  neutron  signal.  As 
discussed  previously  in  Sec.  4,  the  bremsstrahlung  spectrum  incident  on  the  DU  plate  is 
calculated  with  LSP/ITS,  and  MCNPX  is  used  to  calculate  the  neutron  incident  on  the  scintillator. 
This  MCNPX  calculation  is  divided  into  two  steps  to  reduce  computation  time.  In  the  first  step, 
irradiation  of  the  DU  by  the  computed  bremsstrahlung  spectrum  is  used  to  compute  the  output 
neutron  spectrum  from  photofission  and  (y,n)  reactions.  The  second  step  determines  the  time- 
dependent  neutron  energy  spectrum  at  the  end  of  the  hole  in  the  2.43-m  thick  concrete  wall 
(shown  in  Fig.  10).  A  neutron  spectral  resolution  of  100-keV  is  sufficient  to  resolve  the  6.17-m 
TOF  into  25-ns  bins  with  a  statistical  precision  better  than  5%.  The  time  profile  of  neutron 
emission  at  the  DU  is  assumed  to  be  a  trapezoid  with  a  base  width  of  40  ns:  12-ns  rise,  13-ns  flat 
top,  and  15 -ns  fall.  This  trapezoid  approximates  the  LSP/ITS  bremsstrahlung  intensity  variation 
of  photons  above  5  MeV  incident  on  the  DU.  Neutrons  are  tallied  crossing  the  rear  surface  of  the 
lead  shield  at  the  end  of  the  2.43-m  long  hole.  No  attempt  is  made  to  model  the  detector  response 
in  MCNPX.  This  calculation,  in  which  a  bremsstrahlung  spectrum  associated  with  0%  ions  was 
used,  indicates  that  2.35  x  10s  neutrons  are  emitted  into  the  2rc  steradians  centered  about  the 
direction  of  the  15 -cm  diam  hole  in  the  wall  leading  to  the  TOF  detector. 
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Calculated  average 
neutron  spectra  arriving  at 
the  TOF  detector  during 
four  different  25-ns  time 
bins  are  shown  in  Fig.  12. 
The  most  energetic 
neutrons  arrive  at  the 
detector  first,  as  expected. 
However,  all  of  the 
spectra  include  a  cluster 
of  low-energy  neutrons 
that  result  from  scattering 
in  lead  and  concrete 
during  transport.  Only 
about  2700  neutrons 
integrated  over  time  and 
energy  are  transmitted 
through  the  lead  into  the 
vicinity  of  the  neutron 
detector.  About  300 
neutrons  are  incident  on 
the  scintillator  because  the 
scintillator  diameter  is 
only  5.1-cm  compared  with  the  15-cm-diameter  hole  in  the  concrete  wall.  The  number  of 
neutrons  detected  is  even  smaller  because  of  the  efficiency  of  (n,p)  scattering  interactions  in  the 
scintillator.  This  small  number  of  predicted  neutrons  is  consistent  with  the  TOF  detector 
measurement. 

The  calculated  neutron  spectra  averaged  over  25-ns,  Si(En),  exemplified  by  the  spectra  in  Fig. 
12,  are  used  to  compute  the  shape  of  the  TOF  detector  signal  for  comparison  with  a  measurement. 
This  comparison  is  quite  challenging  because  of  to  the  complex  nature  of  the  measured  TOF 
detector  signal.  The  signal  itself  consists  of  piled-up  pulses  that  cannot  be  deconvolved  into 
countable  pulses,  and  even  for  discrete  pulses,  the  proton  recoil  mechanism  and  finite  detector 
size  severely  limit  the  utility  of  using  the  deposited  energy  into  the  scintillator  (pulse  height)  to 
infer  the  incident  neutron  energy  for  that  event.  However,  for  a  detector  operating  strictly  in  the 
linear  regime,  for  fast  scintillation  and  light  collection  times  compared  to  the  signal  timescale, 
and  in  the  limit  of  many  neutron  detection  events  within  the  scintillator  (such  that  an  average 
energy  deposition  can  be  used  to  model  each  neutron  detection  event),  we  make  the  simplifying 
assumption  that  the  detector  signal  can  be  interpreted  as  a  current-mode  measurement.  Assuming 

that  the  detector  response  is  proportional  to  the  neutron  energy,  then  js,(E„)E„dE„  is 

proportional  to  the  number  of  neutrons  detected  in  the  ith  time  interval.  Under  these  assumptions, 
the  current-mode  measurement  describes  the  neutron  energy  deposition  as  a  function  of  time  and 
is  proportional  to  the  computed  energy  deposition. 

All  of  these  weighted  MCNPX  spectrum  integrals  are  combined  to  produce  a  computed 
approximate  time  history  of  the  measured  neutron  TOF  signal.  The  normalized,  computed  time 
history  so  obtained  is  compared  with  the  measured  TOF  signal  for  shot  899  in  Fig.  13.  The 
experimental  signal  has  been  shifted  by  20  ns  relative  to  Fig.  1 1  to  account  for  the  x-ray  flight 
time  from  the  DU  plate.  The  5.08  cm  dia.  by  5.08  cm  long  scintillator  satisfies  the  detector  speed 
requirement.  The  detector  linearity,  which  may  ultimately  be  a  reasonable  approximation,  is 
presently  being  investigated  using  neutron  detection  data  taken  from  a  Van  de  Graaff  accelerator. 
The  large  event-number  requirement  within  each  time  window  is  more  difficult  to  assess.  An 
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Figure  12.  MCNPX-calculated  neutron  spectra  for  neutrons  arriving  at 
the  TOF  detector  during  the  indicated  time  intervals.  Time  zero  is  at  the 
beginning  of  prompt  neutron  emission  from  the  DU. 
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Figure  13.  The  calculated  time  history  of  the  TOF  response 
(green)  from  the  weighted  MCNPX-spectrum  integrals  is 
compared  with  the  measured  response  (red)  for  the  DU  plate 
on  shot  899,  under  several  simplifying  assumptions.  The 
calculation  is  normalized  to  the  measurement. 


average  detector  signal  over 
several  shots  could  help  satisfy 
this  condition  if  more  shots  were 
available  to  analyze.  The 
calculated  neutron  signature  is 
consistent  with  the  measured 
neutron  detector  response,  and 
suggests  that  a  current  mode 
approximation  is  not  out  of  the 
question  for  this  diagnostic. 

Results  show  that  a  prompt 
neutron  signature  with  an 
appropriate  TOF  delay  is 
measured  when  the 

bremsstrahlung  source  induces  a 
fission  response  from  the  DU 
target.  Under  a  set  of  simplifying 
assumptions,  the  computed  time 
history  is  consistent  with  the 
measured  TOF  signal.  When  a 
lead  surrogate  replaces  the  DU 
target,  the  fission  signature 
vanishes,  as  expected. 
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6.  DETECTION  STUDIES  WITH  MASS  ATTENUATION 


Experiments  were  carried  out  with  various  thicknesses  and  types  of  material  interspersed 
between  the  DU  plate  and  the  Mercury  bremsstrahlung  source.  The  materials  were  lead,  2% 

borated  polyethylene  (BPE),  and 
a  combination  of  both  lead  and 
BPE  to  determine  the  impact  on 
the  fission  signature  when  the 
attenuation  was  a  result  of  high- 
Z  and  low-Z  materials.  Both 
delayed  y-rays  and  delayed 
neutrons  were  measured  as  a 
function  of  the  material 
thickness  to  account  for  the  fact 
that  lead  attenuates  y-rays  more 
efficiently  than  neutrons,  while 
BPE  attenuates  neutrons  more 
efficiently  than  y-rays.  As 
stated  in  Sec.  2,  the  AK  gap  for 
these  measurements  was  40  cm. 
The  arrangement  of  the 
attenuation  material,  the  DU 
plate,  and  the  y-ray  and  neutron 
detectors  is  presented  in  Fig.  14. 
The  DU  is  located  2.5  m  from 
the  bremsstrahlung  converter 
and  is  aligned  at  45°  so  that  the 
two  radiation  detectors  view  the 
DU  plate  at  near  normal  incidence.  Delayed  neutrons  were  measured  using  a  single  He-3  detector 
(described  in  Sec.  4)  located  2.0  m  from  the  DU.  Delayed  y-rays  were  measured  using  the  two 
large  NaI:TI  detectors  (40.64-cm  x  10.16-cm  in  area  x  5.08-cm  thick)  placed  3.6  m  from  the  DU 
and  filtered  in  front  with  0.63-cm  thick  lead  and  5.08-cm  thick  BPE. 

Measurements  were  made  with  three  different  types  of  attenuating-material 
configurations:  lead  only,  BPE  only,  and  both  lead  and  BPE.  Thicknesses  were  1.27,  3.81,  6.35 
and  10.16  cm  for  lead,  and  5.24,  45.72,  and  76.20  cm  for  BPE.  For  mixed  lead  and  BPE,  the 
thicknesses  were  1.27-cm  Pb  +  15.24-cm  BPE,  3.81 -cm  Pb  +  45.72-cm  BPE,  and  6.35-cm  Pb  + 
76.20-cm  BPE  with  the  lead  always  adjacent  to  the  DU  plate.  These  thicknesses  are  the  actual 
thicknesses  measured  normal  to  the  DU  plate;  they  are  not  the  thickness  traversed  by  the 
interrogating  bremsstrahlung  (see  Fig.  14).  In  all  cases,  the  attenuating  material  was  located  on 
the  front  and  back  of  the  DU  plate,  as  in  a  sandwich,  and  extended  beyond  the  edge  of  the  DU  so 
that  no  bremsstrahlung  is  directly  incident  onto  the  DU.  The  counts  recorded  by  each  detector 
were  summed  over  a  time  interval,  and  this  sum  is  treated  as  a  single  measurement.  For  each 
configuration,  two  shots  were  taken  and  the  results  averaged. 

Delayed  y-ray  and  delayed-neutron  measurements  are  presented  in  Figs.  15  and  16, 
respectively,  as  a  function  of  attenuation-material  thickness  in  units  of  areal  mass  density,  i.e.,  px 
in  g/cm2,  where  p  is  the  material  density  and  x  is  the  thickness  traversed  by  the  incident 
bremsstrahlung  (  V 2  times  larger  than  the  thicknesses  quoted  above  because  of  the  45°  tilt  of  the 
DU  plate  relative  to  the  machine  axis).  No  attenuating  material  is  used  for  the  zero  thickness  data 
in  Figs.  15  and  16.  The  uncertainty  in  each  measurement  is  typically  smaller  than  the  size  of  the 
data  point  (a  vertical  solid  line  is  used  otherwise).  For  the  delayed  y-ray  measurements  in  Fig.  15, 


He-3 

detector 


Figure  14.  Top  view  of  the  bremsstrahlung  diode,  the  Nal 
detectors  and  the  He-3  detectors  for  measurements  with 
attenuating  material  on  the  DU. 


16 


only  y-rays  recorded  from  3  s  to  5  s  after  the  bremsstrahlung  pulse  and  with  energy  between  3  and 
7.5  MeV  are  included.  This  energy  cutoff  is  chosen  to  avoid  the  208T1  y-ray  line  at  2.614-MeV,  as 
described  in  Sec.  3.  The  3-s  start-time  for  summing  the  counts  is  when  the  Nal  detectors  have 
fully  recovered  from  the  initial  bremsstrahlung  pulse.  (Subsequent  work  has  led  to  a  much 
shorter  recovery  time.)  For  the  neutron  measurements  in  Fig.  16,  pulses  are  summed  from  1  ms 
to  1  s  after  the  bremsstrahlung  pulse.  Again,  the  1-ms  start  time  is  set  by  detector  recovery  from 
the  bremsstrahlung  pulse.  Average  active  backgrounds,  measured  with  no  DU,  are  shown  as 
dashed  horizontal  black  lines  in  Figs.  15  and  16.  These  backgrounds  include  eight  shots:  two 
with  only  lead,  two  with  only  BPE,  two  with  only  a  beryllium  plate,  and  two  with  no  object  or 

attenuating  material.  The 
uncertainty  in  the 
backgrounds  is  the  standard 
deviation,  aB,  for  these  eight 
shots  and  is  indicated  by  the 
shaded  regions  in  Figs.  15 
and  16  (around  the  average 
background  dashed  line). 

The  y-ray  results  in  Fig. 
15  for  all  three  attenuation 
configurations  follow  a 
similar  exponential  decrease 
with  increasing  areal  mass 
density.  The  measured 
decrease  in  y-ray  intensity 
with  increasing  areal  mass 
density  is  only  weakly 
dependent  on  the  atomic 
number  of  the  attenuating 
material.  The  neutron  results 
in  Fig.  16  demonstrate  that 
the  low-Z,  BPE  attenuates 
neutrons  more  efficiently 
than  the  high-Z  lead  for  the 
same  areal  mass,  as  expected. 
The  similarity  of  the  results 
with  BPE  only  and  the  mixed 
lead  and  BPE  is  consistent 

with  neutron  attenuation  being  dominated  by  the  BPE. 

The  LSP/ITS  and  MCNPX  codes  were  used  in  a  manner  similar  to  that  described  in  Sec.  4  to 
calculate  delayed  y-ray  and  delayed-neutron  intensities  with  either  lead  or  BPE  sandwiching  the 
DU  plate.  Results  of  these  calculations  are  presented  as  dashed  lines  in  Figs.  15  and  16  and  can 
be  compared  with  the  measured  counts.  The  estimated  ±50%  uncertainty  in  the  MCNPX 
simulations  is  represented  by  vertical  dashed  lines  at  each  px  that  was  simulated.  This 
uncertainty  is  primarily  from  uncertainties  in  the  238U  photofission  cross  section  between  5  to  8 
MeV  (see  Sec.  4).  For  the  y-ray  results  in  Fig.  15,  the  calculated  intensities  are  normalized  to  the 
measured  counts  with  no  attenuating  material  because  the  absolute  response  of  the  Nal  detector 
has  not  been  determined.  It  is  interesting  to  note  that  the  MCNPX  simulations  suggest  a  bigger 
difference  in  attenuation  between  high-  and  low-Z  material  than  the  measurements  suggest.  This 
discrepancy  is  has  not  been  resolved.  For  the  neutron  results  in  Fig  16,  the  He-3-detector 
absolute  response  is  known’9  and  the  MCNPX  results  are  convolved  with  this  response  to  obtain 
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Figure  15.  Delayed  y-ray  measurements  and  results  of  MCNPX 
simulations  for  the  different  attenuating  materials  as  a  function  of 
the  material  thickness  expressed  as  areal  mass  density,  px,  where  x 
is  the  distance  traversed  through  the  DU  by  the  interrogating  beam. 
The  dashed  lines  are  results  from  MCNPX  simulations  (normalized 
to  zero  px).  The  uncertainty  in  the  data  is  within  the  size  of  the 
symbol  except  where  a  solid  vertical  line  appears.  The  dashed 
vertical  lines  represent  the  +50%  confidence  level  for  the 
simulations. 
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the  calculated  data.  The  neutron  counts  determined  in  this  way  for  both  the  lead  and  BPE 
attenuators  are  remarkably  consistent  with  the  measurements.  The  difference  between  the 

calculated  and  measured 
counts  at  zero  thickness  in  Fig. 
16  is  still  being  investigated. 

As  the  thickness  of 
attenuating  material  is 
increased,  the  y-ray  and 
neutron  fission  signatures 
decrease  so  that  detecting  the 
presence  of  DU  becomes  more 
difficult.  For  the  purposes  of 
this  work,  three  standard 
deviations  above  the  active 
background  (3gb),  represented 
by  the  dashed  red  horizontal 
lines  in  Figs.  15  and  16,  is 
used  as  a  simplified  measure 
of  detection  effectiveness  for 
this  idealized,  laboratory 
setup.  More  sophisticated 
statistical  treatments  will  be 
carried  out  in  future  work  for 
more  realistic  active -detection 
scenarios.  In  Fig.  15,  the 
delayed  y-ray  count  exceeds 
3gb  for  px  up  to  100  g/cm2  for 
all  three  types  of  attenuating- 
material  configurations.  In 
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Figure  16.  Delayed  neutron  measurements  and  results  of  MCNPX 
simulations  for  the  different  attenuating  materials  as  a  function  of 
the  material  thickness  expressed  as  areal  mass  density,  px,  where 
x  is  the  distance  traversed  through  the  DU  by  the  interrogating 
beam.  The  dashed  lines  are  absolute  results  from  MCNPX 
simulations.  The  uncertainty  in  the  data  is  within  the  size  of  the 
symbol  except  when  a  solid  vertical  line  appears.  The  dashed 
vertical  lines  represent  the  +50%  confidence  level  for  the 


Fig.  16,  the  delayed  neutron  count  with  lead  exceeds  3gb  for  px  up  to  160  g/cm  .  Neutron 
measurements  with  BPE  or  mixed  attenuating  materials  require  a  more  detailed  analysis  because 
the  sandwich  nature  of  the  attenuating  material  allows  fission  neutrons  emitted  through  the 
uncovered  edge  of  the  DU  plate  to  be  room-scattered  and  reach  the  detector.  In  the  future, 
ongoing  analysis  of  this  case  will  be  augmented  by  experiments  where  attenuating  material 
completely  surrounds  the  DU  plate. 

Results  show  that  attenuation  with  both  high-  and  low-Z  materials  reduces  the  delayed  y-ray 
and  delayed  neutron  signal  produced  from  bremsstrahlung  irradiation  of  a  DU  object.  The 
attenuation  of  the  measured  delayed  y-ray  signal  is  not  very  sensitive  to  the  Z  of  the  material  and 
depends  mostly  on  the  areal  mass  density;  whereas  attenuation  of  the  delayed  neutron  signal  is 
considerably  larger  with  low-Z  material  than  with  high-Z  material.  Simulations  agree  well  with 
the  measured  neutron  results.  However  simulations  suggest  the  attenuation  for  delayed  y-rays  is 
more  sensitive  to  the  Z  of  the  attenuating  material  than  indicated  by  the  measurement,  with  the 
calculated  attenuation  by  the  low-Z  material  being  somewhat  less  than  the  measurements  suggest. 
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7.  DIODE  IMPROVEMENTS 


To  maximize  the  forward-going  bremsstrahlung  that  induces  fission  in  DU,  electrons  should 
impinge  on  the  anode  converter  at  as  close  to  normal  incidence  as  possible  to  take  advantage  of 
the  preferentially  forward  directed  bremsstrahlung  (half  angle  ~  1/y,  where  y  is  the  relativistic 
factor  for  electrons).  Also,  their  deposited  charge  at  greater  than  about  6-MeV  energy  should  be 
as  large  as  possible.  Both  are  of  particular  importance  for  long  standoff  applications,  where  the 
inverse-square -with-distance  penalty  can  be  very  serious.  In  separate  dedicated  Mercury 
experiments,  the  fission  yield  from  the  DU  plate  was  maximized  by  varying  the  diode  AK  gap 
between  23  and  43  cm  and,  for  the  larger  gaps,  by  adding  a  25-cm-long  insert  to  reduce  the  inner 
diameter  of  the  outer-conductor  wall  upstream  of  the  converter.17' 43  As  the  AK  gap  increases  the 
beam  radius  increases  but  the  angle  of  incidence  becomes  closer  to  normal  on  the  converter.  The 
insert  (not  shown  in  Fig.  1,  but  illustrated  in  Fig,  17c)  reduces  the  beam  diameter  at  the  converter 
while  preserving  near-normal  incidence.17'43  This  insert  was  first  described  and  used  in  Ref.  44 
and  was  called  an  indented  anode  therein.  An  identical  suite  of  diagnostics  is  fielded  for  each 
diode  configuration.  New  current  monitors  on  the  anode  at  the  downstream  end  of  the  insert  (just 
before  the  anode  end  plate)  determine  the  current  actually  reaching  the  bremsstrahlung  converter. 
TLDs  and  an  x-ray  pinhole  camera  are  used  to  measure  the  angular  and  radial  x-ray  dose 
distributions,43,45  and  He-3  detectors  are  used  to  measure  fission  yields  induced  in  the  DU  plate 

both  on-axis  and  15°  off- 
axis.43  All  of  these 
measurements  are  in  broad 
agreement  with  LSP/ITS 
simulations  of  the  electron 
dynamics  and  radiation 
transport  for  each  diode 
configuration,  which 

provides  fundamental 
understanding  of  the 
processes  contributing  to 
fission-yield  variations.43 

Simulations  of  the 
electron  flow  in  the 
Mercury  diode  are 
illustrated  in  Fig.  17. 
There,  LSP-computed  rB0 
contours  (i.e.,  current  flow 
contours  for  azimuthally 
symmetric  current  flow)  of 
the  electron  flow  are 
shown  for:  (a)  a  23-cm  gap 
without  the  insert,  (b)  a  40- 
cm  gap  without  the  insert,  and  (c)  a  40-cm  gap  with  the  insert.  A  stimulated  ion  emission  fraction 
of  0.4%  was  used  in  all  these  simulations. 
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Figure  17.  LSP  rBe  contours  showing  the  enclosed  current  associated 
with  electron  flow  in  the  Mercury  diode  for  a)  a  23-cm  gap  without  an 
insert,  b)  a  40-cm  gap  without  an  insert,  c)  a  40-cm  gap  with  an  insert. 
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In  Fig.  18,  the  measured  current  entering  the  diode  region  (i.e.,  Iload  of  Fig.  2)  is  compared 

with  the  measured  current  actually 
reaching  the  bremsstrahlung  converter  for 
35-,  40-,  and  43-cm  AK  gaps.  The 
plotted  converter  currents  are  four-shot 
averages  (typical  variation  is  less  than 
+10%)  for  each  gap  and  all  shots  have  the 
anode  insert  installed.  The  plotted  total 
current,  shown  for  a  35-cm  gap  (shot 
1313),  is  nearly  identical  to  the  total 
currents  for  the  other  gaps.  LSP 
simulations  of  the  time -dependent 
electron  flow  in  the  diode  for  the  different 
gaps  are  in  agreement  with  the 
experimental  results  in  Fig.  18.  The 
delay  in  current  rise  at  the  converter 
compared  to  the  total  diode  current  is  a 
result  of  early  time  electron  losses  to  the 
outer  conductor  wall  and  to  the  insert  as 
the  e-beam  transitions  from  an  initially 
radial  flow  (to  the  outer  wall)  to  an  axial 
flow  (to  the  bremsstrahlung  converter  and 
anode  end  plate).  At  peak  current, 
however,  the  e-beam  impinges  on  the 
anode  at  larger  radius  with  increasing 
gap,  so  that  some  current  loss  to  the  outer 


Figure  19.  Electron  current  on  the  converter  for  35-, 
40-,  and  43-cm  AK  gaps  compared  to  the  total  current 
entering  the  diode  region.  For  each  gap,  the  current  is 
a  four-shot  average  and  all  shots  have  the  indent 
installed.  The  total  current  is  essentially  the  same  for 
all  gaps. 

wall  and  insert  persist  for  AK  gaps 
larger  than  35  cm,  resulting  in 
lower  peak  current  on  the 
converter  for  larger  gaps.43  This 
explanation  is  supported  by  x-ray 
pinhole  images  which  show  that 
the  radial  distribution  of 
bremsstrahlung  from  the  converter 
shifts  to  larger  radius  as  the  gap 
increases.43 

X-ray  dose  measurements  also 
support  thrs  interpretation.  ’ 

Measurements  of  angular  dose 
distributions  at  1  m  from  the 
converter  show  a  narrowing  of 
angular  spread  with  increasing  AK 
gap,  and  a  corresponding  increase 
in  on-axis  dose  as  illustrated  in 
Fig.  19.  The  insert  is  installed  on 
all  shots  in  Fig.  19  except  the  23- 
cm-gap  shots.  For  each  gap,  an 
average  is  taken  over  4  or  5 
sequential  shots  and  displayed  in 
the  plot  starting  with  shot  1295  for 


Figure  18.  Shot-averaged  angular  dose  distributions  for 
various  gaps.  For  each  gap,  the  average  is  taken  over  4  or  5 
sequential  shots.  The  23-cm  AK  gap  data  was  taken  with  no 
insert,  while  all  other  data  were  taken  with  an  insert. 


the  40-cm-gap  case,  shot  1301  for  the  35-cm-gap  case,  shot  1323  for  the  43-cm-gap  case,  and 
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Figure  20.  He-3  counts  in  100  s  for  a  detector  1  m  from  the  DU 
plate.  The  black  line  is  from  LSP/ITS  predictions.  The  23- and 
33-cm  AK  gap  data  and  the  one  38-cm-gap  shot  displayed  as  a 
solid  green  square  were  taken  with  no  insert,  while  all  other 
data  were  taken  with  an  insert.  The  prediction  assumes  8.9xl04 
fissions  per  He-3  count  as  determined  from  the  computations 
described  in  the  text 


shot  1330  for  the  23-cm-gap 
case.  The  behavior  described 
above  is  attributed  to  electrons 
impacting  the  converter  at  angles 
closer  to  normal  as  the  beam 
flows  at  larger  radius.43  This 
trend  continues  until  the  gap 
reaches  43  cm.  At  this  value, 
electron  losses  to  the  outer  wall 
and  insert  are  large  enough  to 
reduce  the  on-axis  dose.  Shot-to- 
shot  variations  of  the  on-axis 
dose  for  the  23-  through  40-cm- 
gap  shots  are  about  5%,  while 
that  for  the  43-cm  gap  is  about 
14%,  demonstrating  decreased 
reproducibility  due  to  variable 
electron-beam  losses  to  the  insert 
and  wall.  The  23-cm-gap  shots 
without  the  insert  show  the 
widest  angular  distribution. 

Variation  in  delayed-neutron 
fission  yield  with  increasing  AK 
gap  is  shown  in  Fig.  20  for  the 
He-3  detector  located  1  m  in 
front  of  the  DU  plate  in  the 
configuration  of  Fig.  5  (i.e.,  to  the  left  of  the  DU  in  the  figure).  The  insert  is  not  used  on  the  shots 
with  23-  and  33-cm  AK  gaps,  but  is  used  on  the  open-square  data  with  gaps  of  35  cm  or  larger. 
There  is  one  3 8 -cm  gap  shot  with  no  insert  (solid-green  square).  The  black  line  in  Fig.  20 
represents  absolute-value  predictions  from  LSP/ITS  computations  for  the  various  diode 
configurations,  assuming  a  0.4%  ion  fraction,  using  the  procedures  outlined  in  the  Sec.  4  to 
convert  fissions  in  the  DU  to  counts  in  the  detector  (i.e.,  each  He-3  count  corresponds  to  8.9xl04 
fissions  induced  in  the  DU  plate).  The  green-square  datum  at  38-cm  AK  gap  (without  an  insert) 
indicates  that  for  larger  gaps  where  the  neutron  yield  is  larger,  the  yield  is  reduced  without  the 
insert.  Increasing  the  AK  gap  from  23  cm  to  40  cm  and  using  the  insert  increases  the  neutron 
yield  by  a  factor  of  three. 

Predictions  of  electron  angular  distributions  impacting  the  converter  are  shown  in  Fig.  2 1  for 
the  LSP/ITS  computations  in  Fig.  20.43  These  distributions  produce  x-ray  angular  divergences 
consistent  with  the  dose  measurements  of  Fig.  19  and,  as  shown  in  Ref.  43,  neutron-yield 
measurements  from  an  off-axis  DU  plate.  Comparison  of  the  dashed  and  solid  35-cm  AK  gap 
distributions  indicates  the  improvement  expected  with  the  insert.  With  increasing  AK  gap,  the 
computed  electron  angular  distributions  narrow  and  approach  the  converter  at  more  normal 
incidence,  and  both  the  on-axis  dose  and  the  DU  fission-neutron  yield  increase.  For  40-cm  gaps 
or  larger,  LSP  electron-impact  angles  on  the  converter  are  within  5°  of  normal.  For  8-MeV 
electrons,  the  x-ray  divergence  does  not  narrow  significantly  for  smaller  electron  angles. 
Electron  losses  to  the  insert  and  wall  continue  to  increase  for  larger  gaps,  however,  explaining  the 
reduced  neutron  yield  at  43  cm  in  Fig.  20.  If  the  current  lost  to  the  insert  and  walls  could  reach 
the  converter,  LSP/ITS  predicts  a  fission  yield  1. 75-times  larger  than  that  achieved  with  a  40-cm 
gap.  Though  wall  losses  with  a  23-cm  gap  are  small,  the  predicted  20°  electron  impact  angles 
produce  only  about  1/3  of  the  on-axis  fissions  of  the  40-cm  gap.  Therefore,  research  has  begun  to 
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Figure  21.  LSP  predictions  of  electron  angular  distributions 
impacting  the  converter  for  various  AK  gaps  and  diode-load 
configurations. 


modify  smaller-gap  diodes  to 
magnetically  steer  electrons  to 
converter  impact  angles  closer  to  the 
normal.46 

Results  show  that  the  beam 
radius  increases  and  the  electron 
angular  distribution  becomes  more 
normal  to  the  converter  as  the  AK 
gap  increases.  For  a  given  AK  gap, 
the  insert  also  decreases  the  beam 
radius  and  narrows  the  angular 
distribution  moving  it  closer  to 
normal  incidence.  However,  losses 
to  the  wall  and  insert  eventual  limit 
the  benefit  of  these  variations.  By 
increasing  the  AK  gap  from  the 
original  23  cm  to  40  cm  and 
installing  the  insert  described  above, 
an  increase  of  a  factor-of-three  is 
achieved  in  the  bremsstrahlung- 
induced  fission  yield  from  the  DU 
plate. 
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8.  SUMMARY  AND  CONCLUSIONS 


The  IPAD  approach  has  several  features  that  make  it  attractive  for  application  to  the  active 
detection  of  fissile  material.  In  the  work  reported  here,  we  have  demonstrated  that  a  short, 
intense,  bremsstrahlung  pulse  delivered  to  a  DU  plate  at  moderate  endpoint  energy  can  result  in 
favorable  conditions  for  measuring  a  variety  of  bremsstrahlung-induced  fission  signatures  with 
acceptable  active  background  and  low-to-no  passive  background,  even  through  a  moderate-to- 
substantial  amount  of  low-Z  and  high-Z  attenuating  material.  We  present  early  experimental 
results  using  the  8-MeV,  200  kA,  50-ns  Mercury  IVA  at  NRL,  as  well  as  results  from  analyses  of 
these  experiments  and  comparisons  with  simulations.  The  measurements  and  simulations  are  in 
general  agreement,  especially  when  more  recent  photo-fission  cross  sections  are  used.  These 
measurements,  simulations,  and  analyses  demonstrate  the  principles  of  the  IPAD  approach. 
Demonstration  of  principles  is  not,  however,  demonstration  of  a  capability.  Real  world 
conditions  must  be  accounted  for  in  any  real  application.  In  particular,  the  ~  20  rad  (CaF2)  dose 
could  be  considered  prohibitively  high  for  some  applications.  In  that  case,  further  work  would  be 
required  to  take  full  advantage  of  this  IPAD  approach  that  has  the  potential  for  the  least  amount 
of  dose  for  given  detection  scenarios  (be  it  for  short-range  or  long-range  applications)  and  given 
confidence  levels. 

Three  different  fission  signatures  from  a  44-kg  DU  plate  have  been  unambiguously 
measured:  prompt  neutrons  using  a  fast  TOF  detector,  delayed  neutrons  using  appropriately 
filtered  He-3  detectors,  and  delayed  y-rays  using  BGO  and  Nal  detectors.  Additionally, 
measurements  of  delayed  neutrons  and  y-rays  have  been  carried  out  with  three  types  of 
interspersed  attenuating  material:  low-Z  (BPE),  high-Z  (lead)  and  mixed-Z  (BPE  and  lead).  For 
our  idealized,  laboratory  setup,  the  measured  delayed-y-ray  counts  equal  or  exceed  3  times  the 
background  standard  deviation  (3gb)  for  px  up  to  100  g/cm2  for  all  three  types  of  attenuating- 
material  configurations.  The  delayed  neutron  counts  with  lead  equal  or  exceed  3gb  for  px  up  to 
160  g/cm2.  Analyses  and  simulations  of  the  measurements  are  in  general  agreement  with  these 
experimental  results  as  well.  Work  is  underway  to  improve  on  these  results  by  decreasing  the 
detector  recovery  time,  thus  being  able  to  access  any  earlier  and  stronger  fission  signature  that  is 
not  dominated  by  the  induced  background. 

Finally,  we  have  increased  the  fission  yield  by  a  factor  of  3  with  no  change  to  the  driving 
pulsed  power  by  better  control  of  the  election  angles  of  incidence  on  the  converter. 

Future  work  will  focus  on  varying  the  bremsstrahlung  endpoint  energy  to  explore  the  tradeoff 
between  higher  induced  background  and  higher  fission  signal,  carrying  out  experiments  at  large 
standoffs  (10’s  of  m),47  reducing  the  recovery  time  of  detectors  to  access  more  of  the  higher 
intensity  portion  of  the  fission  signature,  developing  new  methods  for  measuring  prompt 
neutrons,47  developing  new  approaches  for  manipulating  the  e-beam  to  maximize  induced  fissions 
without  increasing  the  e-beam  diameter,46  investigating  pulsed-power  architectures  that  can 
reduce  the  size  and  weight  of  the  pulsed-power  driver,20  and  developing  statistical  approaches  to 
quantify  the  probability  of  detection. 
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